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SUMMARY 


There are three mechanisms by which a vertical magnetic field 
might be generated by a nuclear explosion. Asymmetry of the nuclear 
device might produce a non—vertical dipole moment whose time Variation 
would lead to an electromagnetic field with a vertical magnetic component. 
The magnitude of the vertical component due to this mechanism would 
depend sensitively on device design and orientation at the time of the 
detonation. 

Alternatively, Compton currents produced by explosion gamma rays 
might interact with the earth's field and generate a "back" field tending 
to cancel the earth's field in a region about the explosion. Such a 
canceliation would be equivalent to radiation from a Magnetic dipole and 
therefore could lead to a vertical magnetic field. 

Finally, on a much longer time scale, the blast wave motion of the 
heated (conducting) air surrounding the explosion might exclude the earth's 
field and thus generate a low frequency variation in the vertical component 
of the magnetic field. 

In this report only the latter two mechanisms are treated. It is 
Shown that due to air conductivity and the fact that gamma rays travel 
with light speed, the Compton currents do not exclude the earth's field. 
The blast wave motion does exclude the field and the magnitude of the 


resultant electromagnetic field is estimated. 
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THE ELECTROMAGNETIC SIGNAL DUE TO 
THE EXCLUSION OF THE EARTH'S MAGNETIC FIRLD 
BY NUCLEAR EXPLOSIONS 


There are two possible mechanisms by which the earth's magnetic 
field might be excluded from.a large region about a nuclear explosion. 
First, the Compton-recoil electrons generated by explosion gamma rays 
interact with the earth's magnetic field, producing a "back"' field. Approx— 
imate calculations described here indicate that this mechanism does not 
exclude the field. The other possibility is that the field is excluded due 
to the blast wave motion of heated (conducting) air. Calculations indicate 


that this mechanism is effective and its magnitude is crudely estimated. 


HOW THE GAMMA RAYS TEND TO EXCLUDE THE FIELD 

The "prompt" y rays emitted by a nuclear explosion generate 
Compton-recoil electrons in the surrounding air. The recoil electrons 
move roughly a few meters at sea level, predominantly in an outward 
radial direction. As the electrons move, they are turned by the earth's 
magnetic field, thereby producing an azimuthal current which tends to 


cancel the field in the region near the explosion. 


MECHANISM INHIBITING EXCLUSION 
ee ET ANU DLUN 

The slowing—down of the Compton electrons by the air produces 
ionization and makes the air conducting. The conductivity of the air 


leads to a conduction current which opposes the direct Compton current. 
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The question of whether the conduction current cancels the Compton current 
is complicated because the gamma pulse moves with the speed of light, 
which means that the displacement current cannot be ignored and Maxwell's 


equations must be treated in detail — which is done in the next section. 


THEORY 
Denoting the conductivity by oO and the Compton current by rs 


Maxwell's equations become 


> 1090E 47, >. > 
VxB= = a + —(cE +j), (1) 
+ 1 0B 


g at adistance r (kilometers) from the burst point and atime t after the 


explosion is 


2 
o(t, r) = —— nt, r), (3) 

c 

1 
where approximately 
—r /r t 1 
n(t, r) % 4x 107° —- ny [ at f(t! _r/eje hit a (4) 
. rr Q 


In Eqs. (3) and (4), m is the electron mass, vo is the collision frequency 
of thermal electrons with air molecules, \ is the mean free path (kilometers) 


of the explosion gammarays, Y is the explosion yield in kilotons, 7 is 


1 Multiple scattering of the gamma rays is neglected. 
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D the fraction of Y appearing as gamma rays, 8 is the electron attachment rate 


to O07 , and f(t) is the temporal variation of the gamma-ray emission 


oO 
(S dt f(t) = 1). 


The Compton current rf is 


. —r/d t r—(t—t')v,. (t-t')v/a_, 
j(t, r) % —~2x10"e < 5 ny [ R at" #(t" — , \e 0 v(t—t'), (5) 


rr t-— 
c 


where R is the range of the Compton electrons and v is their velocity. 
Assuming that the electrons are scattered radially away from the burst 


point, 


e (6) 


where the coordinate system has its origin at the explosion point with its 
polar axis along the earth's magnetic field and we consider only those 


. . 2 . 
electrons moving in the plane 0 = 7/2. v, 1s the speed of the Compton 


0 
mc v 


eB 


electrons and p = y is the Larmor radius, B being the local 
magnetic field. 

To obtain a precise solution of Eqs. (3) and (4) involves a prohibitive 
amount of calculation. Wecan, however, obtain a reasonably good 
approximate solution. For this purpose, we neglect effects arising from 
the spherical divergence of the Compton electrons and treat the local 


interaction of the electrons and the magnetic field as if the explosion were 


planar. This is justifiable as long as all spatial changes are confined 


“The "back" field from electrons in other directions is less than that 
from electrons in the plane 9 = w/2. 
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within a region whose radial dimension is small compared to the distance 
from the burst point. The results will be seen to justify this approximation. 


Within this approximation, 


— . —> 1 — 7 

E B ey Fe. (7) 

B-B, = b= bes, (8) 
where B, is the earth's magnetic field. Substituting Eqs. (7) and (8) into 


Egs. (1) and (2) and using the planar approximation, 


9b 1 9E. 4a | 

er cat ec (OH + jy)» 8) 
aE 1 ab 

an SOE (10) 


From Eqs. (3) —(6), we find that o and Jy have the form 
o(t, r) = G(r) g(t—r/c), (11) 
J, {ts r) = G(r) h(t —r/c), (12) 
where G(r) varies slowly with r and may be treated as a constant. 
STEADY STATE SOLUTION. 


In view of Eqs. (11) and (12), b and E may be assumed to depend 


only on t = t—r/e. Eqs. (9) and (10) then become 


_ apt) = AT) 4a[o(r) E(t) + i,(7)], (13) 
dE(+t)  db(t) 
aor ie a 4) 


. 3The equation for F can be disregarded since it does not affect the 
B-field. 


and 
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Solving Eqs. (13) and (14) leads straightforwardly to 
jy (7) 
E(t) = — OG) (15) 
jy (7) 
b(t) = “olrT) + A, (16) 


where A is an arbitrary constant. Causality requires that for tr < 0, 


E(t) and b(t) be zero. A simple argument shows that 


i) 
am OT) = 0, (17) 


Hence, A =0Q and 


j, (7) 
b(t) = $ 


a(t) * (18) 


Equations (15) and (18) for E(t) and b(t) satisfy Maxwell's equations 


and all boundary eonditions. They are therefore the required solutions. 


From Eqs. (5) and (6), it is seen that j,(™) < 0. Thus b(t) < 0 


$ 


and therefore adds to the earth's magnetic field. We have the surprising 


result that — not only do the Compton electrons not exclude the earth's 


magnetic field, but in fact they lead to an increase in the earth's field in 


the region where y(t) *# 0. This unintuitive result is a consequence in part of 


the fact that the electromagnetic field due to the displacement current 


moves with the gamma-ray pulse. [If the current source had a velocity 


v<< C, 


then the more familiar result that the field from a given current 
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opposes the field which generates it (Lenz's Law) would be satisfied (see 


Appendix I). 


Equation (18) is also a consequence of the assumed steady state for 


Ip 


current produces radiation which moves with the spe 


may therefore contribute to b. 


shown to be small. 


EFFECT OF TURNING ON j, 


In order to estimate the effect of turning on the current, we will 


assume o is constant and J é 


and (10), 


Let 


Then 


A 
b = 


co 


A - 
b = I dre LET br, t). 


— 0 


A A 
@ b ab 2 2a 
-——— + —— + = 
5 4nrc at k ob 


ee 


is turned on at t = 0. 


K-02 


A 
_ 4mick), 


—27o0(t—T) sin (t — 1) ‘ Ke" — (200) 


t 
— anick dte 
0 


— (200)° 


Then from Eqs. (9) 


A 


Jy 


and o which ignores turning on the current. However, turning on the 
ed of light and which 


In the next section this contribution will be 


(19) 


(20) 


(21) 


(22) 
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' @ 
Hence, 


t 00 Oj,(r , 7) 
—_ ~2n ot —7) os: DO a _ 
b(r, t) af dte J a an K{r—r l, e(t—T)), (23) 


where 
| blr, oy = 2B OP 2g 
and 


00 . 2 2 
i(|r—r'l, e(t—t)) =f dk cos k(r—r') sin e(t—T) Vk “-(270/c)" (24) 


0 4/ Ke ~ (2n0/0)* 
+ J,(2nio Vi (t—7)2— (r=r')“/e") u(t—t—|r—r'|/c). 


Alternatively, 


t 
9 _ 
b(r, t) = —2n far e Tot PY ig(etet-er, T) (r—ct+ecT, T) 
8) 


pte{t-7) 9 / 2 2, 2 
: -f dr jer , 7) s—7 Jy (2nie (t—1T) -—(r—r')"/e )}. 


r—c(t—T) 
(25) 
Assuming again that 
jy(, t) = jp g(t — r/c) for t> 0, (26) 


co) 
where J dt g(t) = 1 and g(t) = 0 for t < 0, 
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we find 


j, g(t—r/c) t 
0 —2 ot — — 
b(r, t) = ———-———- (l--e ed ) — 27 j fat e an o(t—T} | g(2r—r/o~t) 
o 0 6 
rt+c(t—T) 
/ 2 2, 2 
-f[ dr‘ g(t —r'fe) 2 (2xi¢ (t-— T) —(r—r') ie ) | ‘ 
— oe cee Ort O 
r—c(t--T) 
(27) 
The first term gives the steady state solution, Eq. (18). The remaining 
terms give the effect of turning on the current Jy’ 
To illustrate the behavior of b we shall approximate g(t) by a delta 
function, 6(t}). Thus 
j, 6(t — r/c) _ _ 
b(r, t) = 0 (1—e enoty a er olt r/c) 
Oo Q 
(ttr/c)/2 
_ —2ra{t—T) 2 ( ; ye 2 *) 
2A io, dte a7 Jo arias (t—t) —(r—cet) /c 
Or (28) 
j,. &(t—r/c) _ 
b(r, t) = ——___ (1 —e em Ot) 


(ttr/c)/2 _ _ 
— 27c Jo — I dre 2m o(t—7) Jo (2ricy (tr)? (nr —e1) 7/0" ) 


0 


(29) 


It is straightforward to show that for o(t + r/c) >> 1 the last term of 


—27 ot 


Eq- (29) is Ofe ) (see Appendix II). Thus 
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jo &(t — r/c) 


b(r, t) = ——-——. + Oe“). (30) 


More generally it can be shown that for an arbitrary positive and rapidly 


decaying pulse g(t}, the last two terms of Eq. (27) are O(e 27 ot) . 


Therefore, for constant 0, we have verified that the effect of turning 
on the current is small. 
On physical grounds this conclusion could have been expected since 


the radiation from turning on the current will be attenuated by 


Jp 


propagating through the conducting region. We note, moreover, from 


Eqs. (23) and (24) and more directly from Eq. (29) that 


oO 
J ar blr, t) = 0 (31) 


— 0 


as expected from conservation of magnetic flux. 
Since no relevant new features are introduced by the spatial and 
temporal variation of the conductivity, corrections to Eq. (18) are 


expected in general to be small. 


ESTIMATE OF b(t) 


Substituting Eqs. (3) — (6) into Eq. (18) gives 


4 ; ; ; 
Provided, of course, that o is large in the region where the 
current is turned on. 
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10 
4 Yo (T—T ')vp/r 
dr'e(rt [r~r']—)e ON (pe Tt) 
viv c 
7 1 c0 T-Ric 
DC) = Baye 7 
1 _T 
2x10 fatter) a BIT T') 
0 
(32) 
To obtain a numerical estimate of b(t), we let f(t) ~ e', Yo «© c, and 
y x 2. Then 
v ct/x 
+ — — 
b(n) = — (By + Ib(x)]) S$ S#B Lae Oh —er/h) (33) 
j 4x 10 (s) l-e 
X 
a for 0< 7+ < R/c and 
a v R/x 
+ 1— — 
b(t) = —(B, + Ib(7) 1) ——+ = ase (34) 
4x 10 (<) l-e 
X 
: —8 , 
for Ric < t. Inatime of the order 1/@ % 10 sec, b(t) increases 
from zero to its steady state value of 
“c a+B R/x 
b(t) = —(B, + [b(t)])———- == |1 -—e™ “(1 —R/a) |. (35) 
0 4 2 
4x10 (<) 
r 
x il 8 
At sea level, v.~ 1.5 x10 per sec, B % 0.75 x10 per sec, and 
RX 3x 107 cm. Thus 
b(t) % — 0.033 (B, + [b(r)|) (36) 


or 


b(t) # —0.033 B (37) 


; 


= 
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That is, from Eq. (8), 


B = —1.033 By ea- (38) 


HYDRODYNAMIC EXCLUSION OF THE EARTH'S FIELD 


Surrounding a nuclear explosion is an intensely ionized region in 
which the air conductivity is sufficiently high to freeze in the magnetic 
field. As a result, when the blast wave causes the air to expand, the 
magnetic field is carried along. At the same time the magnetic field 
which is external to the ionized region is unable to penetrate the ionized 
region and is pushed outward. 


The behavior of the magnetic field is described by 


—» 


= oo = vx(2xB-c 248) (39) 


t c 4na 


and 
V-B= 0. (40) 


The details of the air heating and the hydrodynamic motion determine 
v and o as functions of time and space. These functions are so 
complicated that a precise solution of Eqs. (39) and (40) would involve a 
lengthy numerical calculation. Instead, we shall make an approximation 
which is adequate to provide semi—quantitative results. 

We assume that surrounding the explosion there exists an isothermal 


sphere of radius R(t) within which the air is highly ionized and the 
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relaxation time tT of the magnetic field is quite long, that is, 


tT dR(t) . 4 Rit) GRit) 
Rit) dt “~~” a2 dt 


>> 1. (41) 


Outside the isothermal sphere, the relaxation time is short and the 
conductivity can be neglected. With this simple model, the field at a 


distance r > R(t) from the explosion is 


2 3 
si yg tr id _ (RO) R(t) ine 
Bs Boe, 3 Bo (0) | | tt | sin 0 e, (42) 


2 3 
R(0) R(t) > 
—By[1— (29 | E cos Oe, 


where again the polar axis coincides with the earth's field B 


— 
e . 
OZ 


Numerical results on the radiation and hydrodynamic flow for a 


sea level nuclear explosion give 


' 1/3 1/3 
R(t) & E + 120( 75) y°'* meters (43) 
1/3 
Y 
for 0 < py tl§ < 0.05, where t is in seconds and Y in kilotons. For 
0.05 < t/yil 3. 
| 1/3 
R(t) = R ~ 50 Y meters. (44) 
max 


Equation (42) holds until Eq. (41) breaks down, after which the second 
and third terms of Eq. (42) decay to zero. The time constant for this 


decay 1s, approximately, 
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2 
4no TO , (45) 
Cc 
where o> is the conductivity of the isothermal sphere when R(t) = R nax . 
A rough estimate of a is 
T,™ 10+* sec . (46) 
2 
470 )R max 2/3 
Consequently, nr, is of the order of afew times Y seconds. 


Cc 


Since this relaxation time is so long, other phenomena, such as fireball 
rise and radiative and convective cooling, will determine the actual time 


for the magnetic field to relax. 
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APPENDIX I 
If in Eqs. (4) and (5), f = f(t—r/v), where v << c, then 


a(t, r) = F(r) g(t —r/v) ~ o(t) (I. 1) 


(t, r) = F(r) h(t — r/v) (tT), (I, 2) 


; 4 
" "6 
where tT = t-r/v. Substituting E = E(t) and b = b(t) into Eqs. (9) 


and (10), we find 


db. ¥ dE _ fev : 
dt cart C (OB ¥ 44) (T-3) 
dE _  v db 
ayo dr (I. 4) 
Solving 
2 eT 2 T 
B(r) = SEW) A an; (7 Yexp{ Sele far" a(s")} (I.5)° 
1—(v/c)” —0co t—(v/c)" oT! 
and 
b(t) =-— = E(t) + b(— 0). (I. 6) 


If b(—o) < O and Jy <0, then b(t) —b(—o) > 0; that is, the field 
behind the pulse is decreased in absolute value. Thus the field is 


partially excluded by the current — as expected from Lenz's Law. 


The non-vanishing of E(t) for large T is a consequence of the 
steady-state, one-dimensional character of the problem. 


ome 
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APPENDIX II 


Consider 


(ttr/c)/2 


eon 5,(2nioV (t— 7)*—(r—er)“/e ) (4) 


r= SPS ~ 2), (II. 2) 


Then 


I = 2 {etzle THI de a totttr/c)z 5 j(2rioV t” (r/o) yz}. 


(II. 3) 
For a(ttr/c)>> 1, 


r= a ttre eT) Pa e tolttr/c)z 5 j(2nio Ve (ele) * az ) 


+(ett+r) feb (II. 4) 


But the integral in Eq. (II.4) can be evaluated exactly and leads to 


_ —27 ot 
I = ee + (ct + r) O(e } . (II.5) 


Hence, 


(II. 6) 


